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On . Abstract 

-^^ , A lower bound for the ground state energy of a one particle relativistic 

*vj ' Hamiltonian - sometimes called no-pair operator - is provided. 

o' 

G^ ■ 1 Introduction 



The multi-particle analogon of the operator 

&: B^A^ ^Do ^) A+ (1) 

is often used to describe relativistic effects in atoms and molecules, see, e.g., 
[pi pi. Here Dq is the free Dirac operator Dq — cia. -V + l3mc^ and 



X 



2 V Eip) 

with E(p) = {c^p^+m^c^y^'^ is the projection operator on the electron subspace 
of the free Dirac operator in momentum space. The underlying Hilbert space 
of B is the space of square integrable four spinors ip with A^ip = ip. 

The operator B was introduced for the two particle case by Brown and 
Ravenhall |g] (see also [^ Q) to cure the continuum dissolution of the two 
particle Dirac operator with Coulomb interaction: since the Dirac operator is 
unbounded form below the two particle Dirac operator has the whole real line 
as spectrum, the eigenvalues of the one particle operator "dissolve" . B is called 
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no-pair operator or in the description of meson models (usually with different po- 
tentials) "reduced Salpeter operator" ijll] . Recently, its multi particle analogon 
was used in connection with the "stability of relativistic matter" -question pO| . 
Similar to the Dirac operator with Coulomb potential, the no-pair operator 
has a critical nuclear charge Z^- Hardekopf and Sucher 0, ItI observed that Z^ = 
2/[(7r/2 -I- 2/7r)a] with the fine structure constant a :— e^ /{he) and investigated 
the ground state energy of B numerically. They claimed that, as for the Dirac 
equation, the ground state energy vanishes for Z = Zc^, p. 2025]. Evans et al [^ 
proved that the energy (ip^Bip) is bounded from below by aZ{l/Tr — 7r/4)mc^, 
if the nuclear charge does not exceed the critical charge Zc, otherwise it is 
unbounded. Tix [|| improved this bound to mc^ (1 - aZ - 0.002Z/Zc) > 0.09. 
This last result shows the difficulties to obtain accurate numerically results near 
the critical coupling parameter and the need for sharp bounds. In this short 
note we explore the results from |15| by using some rather trivial numerics to 
obtain a reliable lower bound for the ground state energy for all Z < Zc- 



2 The Lower bound 

A short summary of the method to obtain a lower bound for [ip, Bip) given 
in p, na] is provided. Any normalized four spinor in the electron subspace of 
Dq can in be written in momentum space as 

^^^^J_(iEo + Eip))uip)\ 

^^^' N{p) \ cp-(Tu{p) ) ^ ' 

with the normalized Pauli spinor u, the Pauli-matrices ctj and Eq = -B(O), 
N{p) = [2E{p){E{p) + Eq)]^/"^. In || it was shown that the minimizer of 
{ipjBip) is radially symmetric, i.e., the minimizing u is of the form u{p) = 
a(|p|)(-\/47r|p|)~^(l, 0)* where a is normalized and positive. This yields for the 
energy {ip,B^I;) 
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e{p)a{pfdp-—-l / a{p')k{p',p)a{p)dpdp' (3) 

/o ^^ Jo Jo 

with e{p) — E{p), p ^ \p\, k :— ko + ki and 

fc,(p',p):=.,(p')0,(^ + ^)s,b). (4) 

The functions Qi are Legendre functions of the second kind (see ||I3] for 
the notation) and appear here for the same reason as in the treatment of 
the non-relativistic hydrogen atom in momentum space [0, problem 77] and 
^j (p) = •\/l + {—\y/e{p). Because of scaling we have assumed h — m = c = 1 
in (0). Using the Schwarz inequality to introduce trial functions hj{p) > in 



the potential energy part of (-0, Bip) 
a{p')kj{p',p)a{p)dpdp' 
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a{p')sj{p')y/Qj y TTvy V M^y ^^ a{p)sj{p)dpdp' 

is obtained. Note that the operator B is similar to the Herbst operator 
H = ^Jcp-p^ + rn^c^ — e^Z/\x\. The ground state energy of H is obtained by 
minimizing the quadratic form (0) with k = ko, without the square roots and 
with 2Z instead of Z . Since no error is introduced by using the Schwarz inequal- 
ity when hj/sj is proportional to the ground state wave function of (0), this 
suggests a similar trial function than the one that was used by Raynal et al |Q 
for the Herbst operator. They obtained excellent upper and lower bounds for 
the ground state of H. Choosing the trial functions hj essentially as the Fourier 
transform of |a;|''~^e~^l'^l, the lower bound 

iiP,Bij)> sup inf£[Z,/3,/i](p). (5) 

-1</3<1 P>0 

was obtained in llSfl where 



aZ ^/p^ + J^ sin[^arctan(p/^)] / 1 



P sin[(/3 + 1) arctan(p//x)] \ e{p) 



2 V- eip)J P + 2 p-_^3/2^(^(p2+^2)-l/2) 



(6) 



with the associated Legendre functions Pv (see M, p. 1060]). To cancel the 
large momentum terms in (ph /3 > is chosen to fulfill 



2 tan(/37r/2) (i 



aZ 13 (1 -/32)tan(7r/3/2)' 



(7) 



The resulting function £[Z,f3,fi\ is numerically evaluated. The values for the 
parameters (3 and fi for some values of aZ can be found in table || and a plot 
of the energy over the coupling constant is shown in figure H. 
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0.9060 


0.106 


0.9843 





0.9003 


0.215 


0.9526 


0.1 


0.8527 


0.423 


0.8652 


0.3 


0.7500 


0.621 


0.7395 


0.5 


0.5709 


0.808 


0.5552 


0.7 



Table 1: Lower bound for the energy E — infp £\Z, /3, /i] (p) and the correspond- 
ing parameters 
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Figure 1: Lower bound for the energy (t/j, B'\\)) of the no-pair Hamiltonian 
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